home *** CD-ROM | disk | FTP | other *** search
/ CU Amiga Super CD-ROM 19 / CU Amiga Magazine's Super CD-ROM 19 (1998)(EMAP Images)(GB)[!][issue 1998-02].iso / CUCD / Programming / LEDA / source / src / plane / _line.c < prev    next >
Encoding:
C/C++ Source or Header  |  1994-11-16  |  2.2 KB  |  111 lines
  1. /*******************************************************************************
  2. +
  3. +  LEDA  3.1c
  4. +
  5. +
  6. +  _line.c
  7. +
  8. +
  9. +  Copyright (c) 1994  by  Max-Planck-Institut fuer Informatik
  10. +  Im Stadtwald, 6600 Saarbruecken, FRG     
  11. +  All rights reserved.
  12. + 
  13. *******************************************************************************/
  14.  
  15.  
  16.  
  17. #include <LEDA/line.h>
  18. #include <math.h>
  19.  
  20. static const double eps = 1e-10;
  21.  
  22. //------------------------------------------------------------------------------
  23. // lines 
  24. //------------------------------------------------------------------------------
  25.  
  26.  
  27.  
  28. line::line()
  29. { PTR = new line_rep; }
  30.  
  31. line::line(const segment& s) 
  32. { PTR = new line_rep(s); }
  33.  
  34. line::line(const point& x, const point& y)    
  35. { PTR = new line_rep(segment(x,y)); }
  36.  
  37. line::line(const point& p, double alpha) 
  38. { PTR = new line_rep(segment(p,alpha,1)); }
  39.   
  40.  
  41.  
  42. bool line::contains(const point& p)   const
  43. { if (vertical())
  44.       return p.xcoord() == ptr()->seg.xcoord1(); 
  45.   else
  46.       return p.ycoord() == y_proj(p.xcoord());
  47.  }
  48.  
  49. bool line::contains(const segment& s) const
  50. { point p = s.start(); 
  51.   point q = s.end();
  52.   if (contains(p))
  53.     if (contains(q)) return true; 
  54.   return false;
  55.  }
  56.  
  57.  
  58. ostream& operator<<(ostream& out, const line& l) 
  59. { return out << l.seg(); }
  60.  
  61. istream& operator>>(istream& in, line& l)  
  62. { segment s; 
  63.   in >> s; 
  64.   l = line(s); 
  65.   return in; 
  66.  }
  67.  
  68.  
  69. bool line::intersection(const line& s, point& inter) const
  70. { double cx,cy;
  71.  
  72.   if (slope() == s.slope()) return false;
  73.  
  74.   if (vertical())
  75.      cx = ptr()->seg.xcoord1();
  76.   else
  77.      if (s.vertical())
  78.         cx = s.ptr()->seg.xcoord1();
  79.      else
  80.         cx = (s.y_abs()-y_abs())/(slope()-s.slope());
  81.  
  82.   if (vertical())
  83.      cy = slope() * cx + y_abs();
  84.   else
  85.      cy = s.slope() * cx + s.y_abs();
  86.  
  87.   inter = point(cx,cy);
  88.  
  89.   return true;
  90. }
  91.  
  92.  
  93. bool line::intersection(const segment& s, point& inter) const
  94. { if (intersection(line(s),inter))
  95.   { double d1  = inter.distance(s.ptr()->start);
  96.     double d2  = inter.distance(s.ptr()->end);
  97.     double l   = s.length();
  98.     if ( d1<=l && d2 <=l) return true;
  99.    }
  100.    return false;
  101. }
  102.  
  103.  
  104. segment line::perpendicular(const point& q) const
  105. { segment s = ptr()->seg;
  106.   point r = q.translate(s.angle()+M_PI_2,-distance(q));
  107.   return segment(q,r);
  108.  }
  109.  
  110.  
  111.